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Automated  System  for  Resolving  the  Position  of  Solar  Radio  Bursts 


'a  swept  frequency  interferometer  system  such  as  the  one  in  operation  at 
Sagamore  Hill  observatory,  described  in  AFGL  Tech  report  76-0194  1976.  contains 
in  its  output  data,  precise  information  which  can  determine,  in  a  plane  resolution, 
the  position  of  a  radio  source  on  the  solar  disk  within  1.  3'  arc.  This  information 
is  the  frequency  separation  between  the  interference  fringes  within  the  spectrum 
observed.  A  system  which  could,  automatically,  convert  this  information  to 
degrees  and  minutes  of  arc  is  described  here,  shown  logically  in  Figure  I. 

A  radio  interferometer  is  a  receiver  system  consisting  of  two  antennas  dis¬ 
placed  from  one  another.  The  system  measures  the  phase  difference  of  signals  re¬ 
ceived  from  the  antennas.  Figure  2  is  a  functional  diagram  of  a  swept  frequency 
interferometer:  X  sin  0  +  d  will  represent  alternately,  a  multiple  of  one  or  a  half 
wavelength  of  the  received  signal.  This  causes  the  phases  of  the  two  signals  to 
alternately  add  and  subtract  at  the  point  of  comparison,  the  hybrid  function.  If  the 
receiver  sweeps  in  frequency  at  a  constant  rate  with  respect  to  time,  the  result  at 
the  output  of  the  receiver  will  be  a  sine  wave  whose  amplitude  is  proportional  to  that 
of  a  discrete  source,  it's  frequence  proportional  to  the  sine  of  angle  o.  The  peaks 
of  this  sine  wave  will  be  referred  to  as  fringes.  In  the  system  at  Sagamore  Hill 
observatory,  the  antennas  are  separated  bv  300  m  on  an  east /west  line.  The 
receiver  sweeps  once  per  second  from  25  MHz  to  75  MHz. 

(Received  for  publication  12  .August  1985) 
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Figure  2.  Sweep  Interferometer 

The  angle.  A,  between  the  source  (SUN)  and  the  baseline  (line  between  antennas) 
changes  as  a  linear  function  of  time.  In  this  system  the  measure  of  time  is  derived 
from  a  stable  10  MHz  oscillator  divided  down  to  0.  1896  Hz  or  I  pulse  every 
3.  27  3  sec;  the  time  period  in  which  the  sun  moves  1.  3 183'  of  arc,  1.  3  18'  arc  is 
a  12-bit  binary  derivative  of  180".  The  0.  1896  Hz  or  5.  273-sec  pulses  are  used 
to  increment  a  14-bit  binary  counter  which  counts  up  to  10,800  in  increments  equiva 
lent  to  1.31831  movement  of  the  sun,  totaling  180’  or  12  hours  of  sun  travel  (see 
Figure  3). 
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Kigure  3.  Time  Base  Generator,  10  MHz  Divided  by 
52736  =  0.  18962  Hz 


The  output  of  the  10,  800  counter,  a  binary  linear  representation  of  time  and 
angle,  is  converted  in  the  sine/cosine  converter  to  a  12-bit  binary  representation 
of  the  cosine  of  the  angle  (TIME),  angle  A  shown  in  Figure  4.  When  the  angle 
from  baseline  to  sun  is  90’  (SUN  OVERHEAD)  the  sign  of  the  cosine  changes  to  a 
negative  vaiue.  If  we  ignore  the  sign  we  have  the  sine  of  an  angle  from  sun  to  a 
plane  normal  to  the  baseline.  This  permits  a  12-hour  positively  incremented 
measure  of  time  and  angle  simplifying  the  process  which  otherwise  would  require 
reversing  the  count  at  noon. 

The  cosine/sine  converter.  12-bit  output,  is  used  to  control,  by  presetting,  a 
I'ounter  counting  lycles  of  the  10 -MHz  clock.  For  a  source  in  the  center  of  the  solar 
disk,  the  output  of  the  counter  then  is  a  signal  whose  frequency  is  inversely  propor- 
lionai  to  the  path  difference  between  the  signals  from  the  antennas  and  proportional 
lo  the  frequency  whose  1/2  wavelength  is  a  sub-multiple  of  the  path  difference. 

This  signal  forms  a  time  base  for  detecting  the  points  in  the  receiver  frequency 
sweeps  where  there  would  be  a  fringe  peak  for  a  solar  disk  center  emission. 


V 


4 


( PATH  DIFF)  X  =  SIN  A'  X  300 

SIN  A' = -COS  A  ORX=-COS  AX300 

ADD  I50m«»ers  FOR  TOTAL  PATH  DIFF. 

Figure  4.  Conversion,  Angle  to  Path  Difference 


The  frequencv  of  tlie  receiver  tuning  is  derived  from  a  tracking  generator  which 
tracks  the  receiver  tuning  continuously  and  produces  a  signal  whose  frequencv,  is 
at  every  moment,  the  same  as  the  receiver  tuning  as  it  sweeps  from  25  to  75  MHz. 

These  two  signals  described  above,  one  whose  frequencv  is  proportioned  to  a 
cosine  function  of  time,  the  other  proportional  to  the  frequencv  observed,  are  com¬ 
pared  in  a  gating  system.  This  system  produces  a  pulse  at  the  exact  point  in  the 
receiver  sweep  where  a  fringe  maximum  would  occur  if  there  were  an  emission 
from  the  center  of  the  sun.  The  difference  between  this  point  and  that,  where  an 
actual  fringe  peak  occurs,  is  proportioned  to  the  displacement  of  the  source  from 
the  center  of  the  solar  disk. 

The  time  base,  a  function  of  the  cosine  of  time  and  also  the  frequency  separation 
)f  fringes  (Af),  will  control  a  gate.  The  gate  can  make  a  measurement  of  the 


frequency  (tracking  generator  output)  every  0.  0004096  sec  or  at  a  rate  of  2441  Hz 
(10  MHz/4096).  However  the  time  it  takes  for  the  sweep  to  change  from  one  fringe 
frequency  position  to  the  next  is  about  0.013  sec  and  the  increment  of  this  change 
to  be  measured  for  resolution  of  1.  3  arc  occurs  in  0.  007  MHz  of  sweep  movement 
in  frequency  domain,  or  at  a  rate  of  35,  7142  Hz  (see  Table  1).  The  problem  in 
making  this  measurement  is  apparent.  One  solution  might  be  to  increase  the 
fundamental  clock  rate  above  10  MHz  to  accomplish  the  desired  resolution,  the 
clock  rate  would  have  to  be  on  the  order  of  7  GHz.  The  state-of-the-art  does  not 
provide  at  this  time  presettable  accurate  digital  dividers  at  much  over  100  MHz. 
A  circuit  which  would  solve  the  resolution  problem  was  conceived  by  never  brought 
to  realization  because  of  time  and  funding  limitations. 


Table  1.  Fringe  Measurements  Limits 

Angle 
of  Source 
Above  Horizon 

1° 

Af 

Separation  of 
Fringes  MHz 

1.  932544 

Af  for 

10° 

MHz 

t  in  Sec 
Sweep  Time 
to  Resolve 

1' 

1' 

A  AF 
Hz 

Number 

of 

Fringes 

10» 

1.484 

0.  448 

1/67K 

746 

33.  6 

20» 

1.  187 

0.  297 

1/lOlK 

495 

42.  1 

300 

1.  000 

0.  186 

1/161K 

310 

50.0 

40- 

0.  875 

0.  125 

1/240K 

208 

57.  1 

500 

0.789 

0.086 

1/349K 

143 

63.3 

600 

0.732 

0.057 

1/526 

95 

68.  3 

700 

0.  694 

0.038 

1/793K 

63 

72.0 

800 

0.  673 

0.  021 

1/1.42CM 

35 

74.  0 

90° 

0.  666 

0.  007 

1/4. 545 

11.  6 

75.  0 

Theoretically  the  circuit  would  make  the  frequency  measurement  of  the  receive: 

sweep  in  the  required  time  with  the  desired  accuracy  by  making 

measurements 

progressively  in 

decreasing  time  and  frequency 

increments  as  the  upward 

moving 

sweep  approaches  each  fringe  position. 

.'Is  shown  in  Figure  5,  the  circuit  would  consist  of  three  reversable  counters. 

One  would  be  counting  cycles  of  the  10  MHz  designated  "T"  counter,  another  would 
count  cycles  of  the  tracking  generator  output  sweeping  from  25  to  7  5  MHz,  designated 
"F"  counter.  The  third  counter  would  also  count  cycles  of  the  tracking  generator 
and  would  be  designated  "AF”  counter. 


Figure  5.  Fringe  Position  Detector 


The  sequence  of  operation  of  the  circuit  shown  in  Figure  5  is  shown  logically 
in  Figure  6. 

The  cosine  divider  divides  the  10 -MHz  clock  frequency  such  that  it's  output  is 
a  series  of  pulses  whose  period  is  proportioned  to  the  cosine  of  time  and  is  exactly 
that  period  in  which  the  cycles  of  the  sweep  frequency  would  fill  the  "F  "  counter 
when  a  fringe  position  is  reached.  The  fringe  positions  occur  in  the  frequency 
domain  at  iiarmonics  of  the  frequency  separation  between  fringes.  The  harmonic 
counter  is  incremented  when  each  fringe  portion  is  detected  and  together  with  its 
coupling  to  the  cosine  counter  maintains  the  pulse  output  period  proportioned  to  the 
occurrence  of  fringe  positions  in  the  receiving  sweeps. 

When  a  fringe  position  is  detected  the  harmonic  counter  is  incremented.  The 
first  pulse  from  the  cosine  divider  following  the  fringe  detection  starts  the  "T" 
counter.  The  next  pulse  from  the  cosine  divider  reverses  the  "T"  counter  and 
starts  the  "F"  counter.  When  the  "T"  counter  reaches  terminal  count  it  reverses 
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itself  and  starts  the  "AF”  counter.  When  the  *  counter  reaches  terminal  count 
it  reverses  itself  and  the  "AF^'  counter  is  held  at  zero  until  the  "t"  counter  reaches 
terminal  count,  then  this  counter  starts  counting.  This  cycle  is  continued  until 
"AF”  reaches  zero  at  the  same  time  as  **T*'  counter  reaches  terminal  count.  At 
this  point  the  fringe  position  has  been  detected.  The  process  continues  through  the 
sweeps.  The  result  is  precise  detection  of  the  fringe  position,  and  the  measurement 
is  made  during  a  period  equal  to  1  cycle  of  the  lO-MHz  clock  or  0.  01  jj.sec. 
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Figure  6.  Logic  of  Fringe  Position  Detector 


The  above  described  method  measures  the  frequency  of  tlic  sweep  and  thi'cefore 
is  independent  of  anv’  non-linearitv  in  the  sweep.  Ollier  nietliods  mav  he  usc*d  to 
^^elect  the  fringe  positions;  some  have  been  tri<*d  using  analog  (’ontrol  or  measure¬ 
ment  of  the  re<’eiving  sweep  but  wer('  n<d  successful  either,  because  of  their  limit^‘d 
accurai’V  or  because  of  the  non-linear  <diaracler  of  most  sw('e[i  rec(d\  er  s. 

'I'he  final  svslcni  output  would  be  a  digital  value  [irofnudioned  to  the  displaretricnt 
)f  the  atdual  solar  emission  from  the  renter  i)f  the  solar  disk  \\’ilh  a  •  or  -  sign  to 
in  iicate  whii'h  -iirection,  east  (»r  west. 
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The  output  circuit  and  logic  is  shown  in  Figures  7  and  8.  An  up/down  counter 
(called  fringe  counter)  would  be  preset  to  its  midpoint  or  1/2  the  value  of  its  terminal 
(full)  count.  The  pulse  derived  from  the  theoretical  fringe  position  for  Solar  disk 
renter  +  90”  would  start  the  counter  upward.  A  pulse  derived  from  detection  of  a 
peak  of  a  fringe  in  the  receiver  system  output  (caused  by  a  solar  emission)  would 
reverse  the  fringe  counter. 
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Figure  7.  Fringe  Read  Out 


The  next  theoretical  fringe  pulse  would  stop  the  count,  read  out  the  number 
and  reset  the  counter  for  the  next  fringe  measurement.  The  number  in  the  counter 
at  the  stop  point  would  be  greater  or  smaller  than  the  count  at  start  (1/2  terminal 
count)  by  an  amount  proportional  to  the  angular  displacement  of  the  real  source  emit¬ 
ting  on  the  surface  of  the  solar  disk  from  the  center  of  the  disk.  Its  sign  (greater 
or  smaller  than  1/2  full  count)  would  indicate  which  direction  east  or  west  that  the 
source  was  displaced  from  the  center  of  the  sun.  This  number  and  a  series  of  these 
numbers  could  be  totalled,  then  averaged.  (Divide  the  number  of  measurements 


, 


to  give  an  integrated  result  of  the  measurements. )  The  extent  of  integration  could 
be  controlled  by  the  number  of  measurements  used  (divider  of  total).  This  would  be 
the  final  output  of  the  system. 
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Figure  8.  Resolver  for  Fringe  Position  Displac  ement 
From  Center  of  Solar  Disk  Logic- 


There  have  been  changes  in  the  state-of-the-art  of  devic-es  which  might  have 
application  in  this  kind  of  task,  perhaps  offering  a  simpler  solution.  Time  and 
funds  were  the  main  factors  in  preventing  final  realization  of  the  conc-ept  described 
above. 

The  same  general  principles  of  the  system  described  could  be  applied  to  a 
resolver  of  Solar  emission  position  using  the  data  from  the  fixed  frequency  measure¬ 
ments  in  the  RSTN  system.  Using  two  antennas  spac-ed  several  wavelengths  (greater 
number  of  wavelengths  spacing  =  greater  resolution)  a  synthetically  introduced 


constant  phase  rotation  of  the  signal  from  one  of  the  antennas  to  a  hybrid  would 
produce  fringes  or  modulation  of  the  receiver  output.  The  frequency  of  modulation 
or  spacing  of  the  fringes  as  they  are  related  to  the  synthetically  introduced  phase 
shift  would  give  the  information  needed  to  resolve  the  source  position.  This  ap¬ 
proach  to  the  problem  might  be  more  fruitful  than  with  the  swept  frequency  data  as 
the  signal  to  noise  ratio  in  the  fixed  frequency  data  is  much  greater. 
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